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if other particle "fog gcn«r-raiinf functions 



TO, invention relates o to fleW of fluorescence s^umopy, and more particularly to a 
Bri* te eing charge **d <"-«-." I"""-" 1 m0 " C1 " eS OT 
panicles present in a sample 

The primary data of an e*penment in fluorescence correlation spectroscopy (FCS) is a 
seq^nce of photon counts delected from a microscopic measurement volume- An essential 
attribute of the fluorescence correlation analysis is the calculation of Uve second order 
autocorrelation function of photon detection. This is a way how a stochastic tuncUon (of 
photon counts) is transformed into a statistic*! function having an expected shape, semi* as a 
me ans to estimate some parameters of the sample. However, fre calculation of the 
autocorrelanon function is not the only way for extracting utformadon about Ote sample from 
the sequence of photon counts. Further approaches are based on moment analysis and analysis 
of the distribution of the number of photon counts per given time interval (Qian and Elson 
Proc. Natl. Acad. Sci. USA. 87: 5479 . 483. 1990; Qian and Elson. Biophys, J. 57: 375 - 380. 

1990). 

4 

*, intensity of fluorescence detecned <ron> ■ ■** » .«■* " * "rt™ * 
«_* on to cognacs of to P^e «oth respect ,0 to focus <* <* ■ 
Wore, a reliable Uq-b- of measurements should account for to |«l «f to 
Ufa** moment volume. Even th«tf> to culculadon of a theoretical tM.< 
O. number of pboton counts is -or. complex k.kH^ 1— *- «" * 
one to MM- * to of photon. —» — ** *»■* « 0f "* 

cnUadon and to M«Mc brightness of fluoreacen, .pedes, * therefore, to 
MM. of to number of photon eonrO. can be «*. for sample ««*. TV J- 
.^ific bn^ess" dene*, to M count - of to de«ctor from W • 

p^icl, of *v.n species situ«ed in . crt* point » to s»ple. c^vendonaU, m to pou« 
where the value of to spatial brightness profile tunc*™ is unity. 

Tie Bra. realign of this kind of analysis fas demonstrated on to oasis of moments of to 
pho»n count number dismbudon (Qian and Bson. Proc. Nad. Acad. Sot. USA. 87: 5*79 - 



483, 1990). The ^ th factorial moment of the photon count number distribution />(*) is 



defined 



as 



a) 



In turn, factorial moments are closely related to factorial cumulants; 



(2) 



or 



' ' /-I 



(3) 



(■<? . are binomial coefficients, and **s are cumulants.) The basic expression used in moment 
analysis, derived for ideal solutions, relate *-th order cumulam to concentrations ( C ,) and 
specific brightness values (q t ) 



(4) 



Here, & is the *-th 



' (V) 




(5) • 



Usually in FCS. the unit of volume and the unit of B are selected which yield % v s Zi = ^ 
After selecting this convention, concentrations in the equations are dimensionless, expressing 
the mean number of particles per measurement volume, and the specific brightness of any 
Species equals the mean count rate from a particle if situated in the focus divided by the 
numeric value of B(0). The value of this constant is a characteristic of optical equipment It 
can be calculated from estimated parameters of the spatial intensity profile (see below) Qian 



and Elson used experimental values of the first three cumulants to determine unknown 
parameters of the sample. The number of cumulants which can be reliably determined from 
experiments is usually three to four. This sets a limit to the applicability of the moment 
analysis. 




The idea behind the so-called fluorescence intensity distribution analysis (FDA), which has in 
detail been described in the international patent application PCT/EP 97/05619 (international 
publication number WO 98/T68 14), can be well understood by imagining an ideal case when a 
measurement volume is uniformly illuminated and when there is almost never more than a 
single particle illuminated at a time, similar to the ideal situation in cell sorters. Undo- these 
circumstances, each time when a particle enters the measurement volume, fluorescence 
intensity jumps to a value corresponding to the brightness of a given type of particles. 
Naturally, the probability that this intensity occurs at an arbitrary time moment equals the 
product of the concentration of a given species and the size of the measurement volume. 
Another fluorescent species which may be present in the sample solution produces intensity 
jumps to another value characteristic of this other species. In summary, the distribution of 
light intensity is in a straightforward way determined by the values of concentration and 
specific brightness of each fluorescent species in the sample solution. 



• 



It is assumed that the light intensity reaching the detector from a particle as a function of 
coordinates of the particle is constant over the whole measurement volume, and zero outside 
it Also, it is assumed that the diffusion of a fluorescent particle is negligible during the 
counting interval T, In this case, the distribution of the number of photon counts emitted by a 
single fluorescent species can be analytically expressed as double Poissonian: the distribution 
of the number of particles of given species within this volume is Poissonian, and the 
conditional probability of the number of detected photons corresponding to a given number of 
particles is also Poissonian. The double Poissonian distribution has two parameters: the mean 
number of particles in the measurement volume, c and the mean number of photons emitted 
by a single panicle per dwell time, qT. The distribution of the number of photon counts n 
corresponding to a single species is expressed as 



(6) 



where m runs over the number of molecules in the measurement volume If PXn) denotes the 
distribution of the number of photon counts from species r, then the resultant distribution P\n) 
is expressed as 



TWs means that /*(/!) can be cadcuJated as a convolution of the series of distributions P\(n). 

■ . ■ ■ ■" >■ ■ ■ ■ - ■)■•'.' 

Like in FCS, the rectangular sample profile is a theoretical model which can hardly be applied 

in experiments. One may divide the measurement volume into a great number of volume 
elements and assume that within each of them, the intensity of a molecule is constant 
Contribution to photon count number distribution from a volume element is therefore double 
Poissonian with parameters cdV and qTB{r). (Here q denotes count rate from a molecule in a 
selected standard position where B = 1, and 5(r) is the spatial brightness profile function of 
coordinates.) The overall distribution of the number of photon counts can be expressed as a 
convolution integral over double Poissonian distributions. Integration is a one-dimensional 
rather than a three-dimensional problem here, because the result of integration does not 
depend on actual positions of volume elements in respect to each other. Figuratively, one may 
rearrange die three-dimensional array of volume elements into a one-dimensional array, for 
example in the decreasing order of the value of B. 

* .. . i ' • T 

In a number of first experiments described in the international patent application PCT7EP 
97/05619, the photon count number distribution was indeed fitted, using the convolution 
technique. The sample model consisted of twenty spatial sections, each characterized by its 
volume Vj and brightness B r -However, the technique described in this patent application is 
slow and inconvenient in cases involving a high number of samples to be analyzed, like in 
diagnostics or drug discovery, or m analyzing distribution functions involving more than a 
single argument. 

* -f r ' ■ - 

■4 ' * ' 
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Therefore, it is an object of the present invention to present a convenient and much faster 
technique for analyzing fluorescence intensity fluctuations. 



According to the present invention there is provided a method for characterizing fluorescent 
molecules or other particles in samples, the method comprising the steps of: 

a) monitoring fluctuating intensity of fluorescence emitted by the molecules or other particles 
in at least one measurement volume of a non-uniform spatial brightness profile by 
measuring numbers of photon counts in primary time intervals by a single or more photon 
detectors, 

b) determining at least one distribution function of numbers of photon counts. P (n), from 
the measured numbers of photon counts, 

c) determining physical quantities characteristic to said particles by fitting the experimentally 
determined distribution function of numbers of photon counts, wherein the fitting 
procedure involves calculation of a theoretical distribution function of the number of 
photon counts P(n) through its generating function, defined as G(|) = />( n ). 



The formal definition of the generating function of a distribution P(n) is as follows: 



m 

C<a = Lf P(n). 



A«0 



(8) 



What makes the generating function attractive in count number distribution analysis is the 
addidvity of its logarithm: logarithms of generating functions of photon count number 
distributions of independent sources, like different volume elements as well as different 
W species, are simply added for the calculation of the generating function of the combined 
distribution because the transformation (8) maps distribution convolutions into the products of 
the corresponding generating functions. 

In a particular preferred embodiment, one might monitor the fluctuating fluorescence intensity 
in consecutive primary time intervals of equal width. Typical primary time intervals have a 
width in the order of several tenth of microseconds. The total data collection time is usually 
several tenth of seconds. 

In a further preferred embodiment, numbers of photon counts { m J subject to determination of 
a distribution function P(n) u» step b) are derived from numbers of photon counts in primary 



time intervals {Nj} by addition of numbers of photon counts from primary time intervals 
according to a predetermined rule. One might e.g. be interested in choosing numbers of 

- V r 

photon counts {n*) subject to determination of a distribution function P(n) which are 
calculated from the numbers of photon counts in primary time intervals {Nj) according to the 

rule n, = £N„ i ^, where M is an integer number expressing how many times the time 
interval in which f «i } is determined is longer than the primary time interval- 

' • ' . . , '.' .'" -'**'" ' - ' .. , .. . ' , • *~ ' ; ' .' .... , > t .. .' " ' .... . .• i ' r .--,i 

In a further embodiment, numbers of photon counts (nj are derived from predetermined 
primary time intervals according to a rule in which primary time intervals are separated by a 

time delay. In particular, me following rule can be applied: n, - £(N„,^ +N^ (J+t)+1 ) r where 

M and L are positive integer numbers, {m} are numbers of photon counts subject to 
determination of a distribution function P(n) t and (Nj) are the numbers of photon counts in 
primary time intervals. 

In some cases, it might be preferred to determine not only a single distribution function in step 
b), but ratiier a set of distribution functions P(n) , These can be determined according to a set 
of different rules, said set of distribution functions being fitted jointly in step c). As an 
example a set of distribution functions with different values of M and/or L might be fitted 
jointly. 

Typical physical quantities which might be determined in step c) avoiding to the present 
invention are concentration, specific brightness and/or diffusion coefficient of molecules or 
other particles. 

" " ■ * ' J" - ■ ■ 1 - " ' " f ' 

In a further preferred embodiment, the generating function is calculated using the expression 

0(4) = exp[ J*?c(<7) J d'ne^'™ - 1)] ; where c(q) is the density of particles with specific 

brightness q, T is the length of the counting interv.al t and B(r) is the spatial brightness profile 
as a function of coordinates. 



Applying the definition (8) to formula (6) with c->cdV and q -* qB{r) , the contribution 
from a particular species and a selected volume clement dV can be written as 



G, (£c/V) = exp[o,<*V(e t » , - l> «' r ' ( " - 1)). 



(9) 



Therefore, the generating function of the total photon count number distribution can be 
expressed in a closed form 



Ot6 = exp(J>J (e't-'W - Y)dV}. 



(10) 



Numeric integration according to Eq. (9) followed by a fast Fourier transform is the most 
effective means of calculating the theoretical distribution P(n) corresponding to a given 
sample (i.e., given concentrations and specific brightness values of fluorescent species). If one 
selects £ = , then the distribution P(n) and its generating Function G(p) are interrelated by 
the Fourier transform. Therefore, it is particularly preferred to select the argument of the 
generating function in the form <f = and to use a fast Fourier transform algorithm in 

calculation of the theoretical distribution of the number of photon counts out of its generating 
function. 

When calculating the theoretical distribution P(n) in step c) according to the present 
invention, the spatial brightness profile might be modeled by a mathematical relationship 
between volume and spatial brightness. In particular^ one might model the spatial brightness 

^^^^ 

profile by the following expression: — = a t x + a t x 2 + a 3 *\ where dV denotes a volume 

element, x denotes logarithm of the relative spatial brightness, and a\ , a 7 and a 3 are empirically 
estimated parameters. 

Some fluorescent species may have a significantly wide distribution of specific brightness. For 
example vesicles, which are likely to have a significantly broad size distribution and a random 
number of receptors, may have trapped a random number of labeled Ugand molecules. In order 
to fit count number distributions for samples containing such kind of species, it is useful to 





m 
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modify Eq: (10) in th e flowing manner. The assumption is made that the disunbution of 
brightness ofparticles/y within a species is mathen^^^ 



(11) 




of convenience: all moments of this distribution 
can be analytically calculated, using the following formula: ^ : 



I 



l_a-M 



(12) 



It is straightforward to derive the modified generating function of" a photon count number 

distribution, .One can rewrite (9) as follows: 




Q(& = expfXe, [dv\dqp{q\a (1 b t Xe 1 *- 0 ?" 1 '? -I)], 



where 



(13) 



"(a) 



(14) 



The integral over s can be performed analytically: 



-1 




(15) 



The parameters a, and A arc related to the mean brightness 4 
distribution of by 



ess 



A f = 



.2 ••; 



(16) 



In the range of obtained count numbers, the probability to obtain a particular count number 
usually varies by many orders 0/ magnitude, see for example the distribution of Fig. 1. 
Consequently, the variance of the number of events with a given count number has a strong 
dependence on the count number. To determine weights for least squares fitting, one may 
assume that, light intensities in all counting intervals are independent. Under this assumption, 
one has a problem of distributing M events over choices of different count numbers n, each 
particular outcome having a given probability of realization. P(n), Covariance matrix elements 

of the distribution can be expressed as follows: 




M 1 



(17) 



where M is the number of counting intervals per experiment. 

For a further simpliBcation, one may ignore the second term on the right side of Eq. (17), 
which can be interpreted as a consequence of normalization. In this case, the weights simply 
equal to the inverse values of the diagonal covariance matrix elements 



W. = 



M 



P(n) 



(18) 




Dispersion matrix (17) corresponds to the multinomial distribution of statistical realizations of 
histograms. The Poissonian distribution, with the constraint that the total number of counting 
intervals M is fixed, will lead to the multinomial distribution. This is the rationale behind 
using Poissonian weights as given in Eq. (18), 



Let n, be the expectation value of the number of events of counting * photons and let 
# = be their sum. Let *j be a statistical realization with M = Z t m t . Assume that 
realizations ir^,/n,,... obey Poissonian statistics 



• v ■ . 

[/FtJ } »/R l t...J: 



(19) 



whbrc ^ have . introduced the . notation n^n* . . .» .. . !] K ~T l and 
m^'.m,\ ... s [/»i 0 , m,,...). The probability of to 



P(M)-——e- s 



The conditional probability of having m^.m^,.:. evehts-if there is a total of Af events is 



or 



This is the multinomial distribution where we have introduced: 



(21) 



and C^ ., " are multinomial coefficients. 

In general, a linear or unearned least squares fitting returns not only the values of the 
estimated parameters, but also their covariewce matrix; provided the weights have been 
meaningfully, set. It may turn out to be possible to express the statistical errors of the estimated 
paruneters analytically in some simple cases (e.g., for the rectangular sample profile and 
single species) but in applications at least two-component analysis is usually of interest 
Therefore, one may be satisfied with the numerical calculations of statistical errors. In 
addition to the ..theoretical" errors with the assumption of non-correlated measurements (Eq. 
(17)), in some cases statistical errors have been estimated empirically, making a series of 
about a hundred FID A experiments on identical conditions. As a rule, empirical errors are 
lugher than theoretical ones by a factor of three to four, Empirical errors appear to be closer to 
the theoretical ones in scanning exp runents. Therefore we are convinced thai the main reason 
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of underestimation of theoretical errors is the assumption of non-correlated measurements. 
Table t compares statistical errors of parameters estimated by fitting a photon count number 
distribution (FIDA) according to the present invention and by the moment analysis. Error 
values are determined through processing a series of simulated distributions, The present 
invention is overwhelmingly better than the moment analysis if the number of estimated 
parameters is higher than three. 



Tablet 







Primary 


Number 


Number of 


" Parameter 


Value 


Percent 


1 Percent 


Data 


time 


of species 


estimated 


specification 


(qs in 


error of 


error of 


collection 


interval. 




parameters 




kHz) 


FIDA 


moment 


time, s 


|ls 










(inventi- 


analysis 














on) 




10.0 


40,0 


1 


2 


c 


0.5 


0.59 


10.54 










9 


60.0 


0.56 


0.51 


10.0 


40.0 


2 


T 




O05 


2.00 


[271 












150.0 


1.54 


1.89 










c ? 


3.0 


0.53 i 


0.62 










m 


5.0 (fixed) 






10.0 


40.0 


2 


4 


Cl 


0.05 


2.26 | 


4.99 










?' 


150.0 


1.63 


2.53 












3.0 


3.18 


17.6 












5.0 


3.35 


14.9 



Confocal techniques may be applied to a wide field of applications, such as biomedicine, 
diagnostics, high through-put drug screening, sorting processes such as sorting of particles 
like beads, vesicles, cells, bacteria, viruses etc.. The conjugate focal (confocal) technique is 
based on using a point source of light sharply focused to a diffraction-limited spot on the 
sample. The emitted light is viewed through a spatial filter (pinhole) that isolates the viewing 



area to that exactly coincident with the illuminating spot. Thus, the illumination and detection 
apertures are optically conjugated with each other. Light originating from focal planes other 
than that of the objective lens is ; rejected, which effectively provides ? very small depuh of 
field. Therefore, in a particular preferred embodiment of the present invention, in step a) a 
confocai microscope is used for monitormg^e intensity of nubres^^^ 
high signal-to-noise ratio, it is useful t P monitor the intensity of fluorescence using an 
apparatus that comprises: a radiation source (12) for providing excitation radiation (14). an 
objective (22) for focussing the excitation radiation (14) into a measurement volume (26). a 
detector (42) for detecting emission radiation (30) that stems from the measurement volume 
(26), and an opaque means (44) positioned in the pathway (32) of the emission radiation (30) 
or excitation radiation (14) for erasing the central part of the emission radiation (30) or 
excitation radiation (14). It might be particularly preferred to use an optical set-up described in 
detail in Figure 9. 

In a further preferred embodiment, the method according to the present invention is applied to 
fit a joint distribution of photon count numbers, to experiments, fluorescence from a 
microscopic volume with a fluctuating number of molecules is monitored using an optical set- 
up (e.g. a confocai microscope) with two detectors. The two detectors may have different 
polarizational or spectral response. In one embodiment, concentrations of fluorescent species 
together with two specific brightness values per each species are determined. The two- 
dimensional fluorescence intensity distribution analysis (2D -FIDA) if used with a polarization 
cube is a tool which can distinguish fluorescent species with different specific polarization 
ratios. This is a typical example of a joint analysis of two physical characteristics of single 
molecules or other particles, granting a significantly improved reliability compared to 
methods focussed on a single physical characteristic. 

In order to express the expected two-dimensional distribution of the number of photon counts, 
it is favourable to use the following assumptions; (A) Coordinates of particles are random and 
independent of each other. (B) Contribution to- fluorescence intensity from a particle can be 
expressed as a product of a specific brightness of the particle and a spatial brightness profile 
function characteristic to the optical equipment. (C) A short counting time interval T is 
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selected, during which brightness of fluorescent particles does not significantly change due to 
diffusion. 





At first, a joint distribution of count numbers from a single fluorescent species and a single 
small open volume clement dVis expressed. The volume element is characterized bv 
coordinates r and spatial brightness B{r), and the fluorescent species is characterized by its 
specific brightness values and q % . By q t and q t , mean photon count rates by two 
detectors from a particle situated at a point where B(r) = 1 are denoted. A convenient choice is 



to 



select a unit of B, as usual in FCS, by the equation Xi = X%. where x t - fV(r)d 3 r. If the 



volume element happens to contain m particles, then the expected photon count numbers per 
time interval T from the volume element are mq t TB(r) and mqJB(r), while the distribution 
of numbers of photon counts from m particles /»(/!„ flj Im) is Poissonian for both detectors 
independently: 



I m) = K ra <^ ^mAmm^^w 



(22) 



From the other side, under assumption (A), the distribution of the number of particles of given 
species in the volume element is Poissonian with mean cdV , c denoting concentration: 



_{cdV) 



m 

iv 



/WW (23) 



The overall distribution of the number of photon counts from the volume clement can be 
expressed using Eqs.22 and 23: 



!* k - 1 ■ * " '** - 1 ■ > * " 

In this formiila, a contribution Erom background count rates; X\ by detector 1 and by 
detector 2; as well as contributions from different fluorescent species, denoted by the subscript 
f; have been integrated. Numeric integration according to Eq. 27 followed by a fast Fourier 
transform is a very efficient means for calculation of the theoretical distribution P^n^n^) 

corresponding to a given sample (i.e. given concentrations and specific brightness values of 
r fluorescent species). 

- ■ ' * - -. ■ ; . - 

The spatial brightness function is accounted through the spatial integration on the right side of 
Eq. 27. The three-dimensional integration can be reduced to a one-dimensional one by 
replacing three-dimensional coordinates r by a one-dimehsional variable, a monotonic 
function of the spatial brightness £(r). A convenient choice of the variable is 

x = ln[£(0) / B(r)] ; A sufficiently flexible model of the spatial brightness profile is presented 

. • •■ > . ' -■ ... , ■ . 

by the following expression: 

■'■ . . - . • ' - ■ ■■;./■(--- 

— -oc^l + fl^+OjX 1 ). ■ . , . , (28) 

■ . ' ■ '/ ■ 

In the interval of obtained count numbers, the probability to obtain a particular pair of count 
numbers usually varies by many orders of magnitude. Consequently, the variance of the 
experimental distribution has also a strong dependence on the count numbers. To determine 
weights for least squares fitting, for simplification it is assumed that coordinates of particles in 

1 * - 

all counting intervals are randomly selected. (This means one ignores correlations of the 



coordinates in consecutive counting intervals.) Under this assumption, one has a problem of 
distributing M events over choices of different pairs of count numbers n lt m, each particular 
outcome having a given probability of realization, P(n lAl ). Covariance matrix elements of the 
distribution can be exprwsed as follows: 



M 



(29) 




where M is the number of counting intervals per experiment. 



For a further simplification, one may ignore the second term on the right side of Eq. (29), 
which can be interpreted as a consequence of nomiaiiation. In this case, the weights simply 
equal to the inverse values of the diagonal covariance matrix elements 



M 



(30) 




In general, a linearized least squares fitting algorithm returns not only values of estimated 
parameters, but also their covariance matrix, provided weights have been meaningfully set. In 
addition to "theoretical" errors corresponding to the assumption of uncorrected measurements 
(Eq. (30)), in some cases statistical errors heve been determined empirically, making a series 
of about 100 2D-FIDA experiments at identical conditions. As a rule, empirical errors are 
higher than theoretical ones by a factor of three to four. The main reason of underestimation of 
theoretical errors is most likely the assumption of uncorrelated measurements. 
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Even though the assumption of uncorrected measurements yields underestimated error 
values, it is a very useful theoretical approximation, allowing 10 compare accuracy of analysis 
under different experimental conditions as well as different methods of analysis. Also, this 
approximation provides an easy method of dau simulation which is a useful tool in general: A 
simple and very fast method of data simulation is calculation of the expected event number as 
a function of photon count numbers and addition of random Poisson noise to each event 
number independently; •• - : • : — • • - 7 ..; — - ■' - . — - ; 



fo Table 2 theoretical errors of one-dimensional and two^dimensibna] fluorescence ihtens^ 

■WY^^pn>are^prcsentW in two selected cases of 
two fluorescent species. In both cases the ratio of spedfic brightness values of the two species 
is three. In the case of 2D-HDX, it is assumed that spectral sensitivities of toe, two detectors 
are tuned to different species. In both cases data collection tune of 10 s, time window of 40 ^ 
and background count rate of 1 kHz are assumed. Note that the statistical errors of the 
estimated parameters are signiflcanUy loWer in the 2D-FTDA example 




.A. 




. 1.* i_ ■ 



mm 



- %*C*£-: 

"■It 



Table 2 



Method Parameter specification Parameter value (selected) Percent error 



1D-FIDA c,V 


0.5 


6.6 


(according cjV 


0.5 


4.9 


to the qi 


60 kHz 


2.2 


invention) q? 


20 kHz 


10.2 


2D-FDA C|V 


0.5 


1.1 


(according CjV 


0.5 


1.1 


to the q A i 


60kHz 


0.T7 


invention) q*j 


20kHz 


1.2 


Obi . 


20 kHz 


1.2 


QB2 


60kHz 


0.77 



The two-dimensional fluorescent intensity distribution analysis according to the present 
invention is in the following compared to the prior art moment analysis. 

Factorial moments of the distribution P{n u m) aze defined as 

Factorial moments are related Lo factorial cumulants Af u 

K u = F u - £ Ct-'C'jK,.,^. (32) 
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C denotes binomial coefficients. Cumulants can be expressed through concentrations and 
specific brightness values by a simple relation 



•(33) 



If the unii of B is selected by the equation^, = has the meaning of the sample volume, 

denoted by V, and Eq. 33 can be written as 



(34) 




where Y denotes a series of constants characterizing the brightness profile: 



(35) 



The principle of moment analysis is to determine valves of a few cumulants from an 
experiment and solve a system of Eqs. 33 in respect to unknown concentrations and brightness 
values.. ■' ;. ' ' ' ■• : v ■ ■:. - ^ ■ ■ 



In Table 3 statistical errors of 2D- FDD A (present invention) and 2D-MAFED (prior art) are 
presented determined by generating a series of 30 random distributions of count numbers, 
simulated for identical "samples", thereafter applying 2D-FIDA and 2D-MAFTD, and 
determining the variance of estimated parameters in both cases. Note a tendency that the 



nnBHi.' 



advantages of 2D-FIDA compared to 2D-MAFID increase with the number of parameters to 



be estimated. 
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Data 


Number 


Number of Specifica- 


Specifi- 


True 


Percent 


Percent 


collec- 


of 


estimated tion of 


cation of 


values 




Vli \Jl XJl 


tion 

%M \f %3 


SDccies 


parameters cumulants 


Ddrame- 




FTDA 




time, s 




used in 
MAFID 


ters 




■ 




2 


1 


J JV/ll 


cV 


05 


1 1 


1 X 






if 


MA 


60 kHz 


1 2 








Km 


43 






1.0 


in 


it 




p< V 




i 1 


2.5 








W2 v 




i n 


4.Z 




■ 


x Ml2 


4A1 




n 7n 
u. /u 










Q A1 


20 kHz 






• 


• 

• • 




Out 


40 kHz 










• 




80 kHz 




\11ACU/ 


10 


3 


5 Kai 


CiV 


0.2 


20 


24 






if 


CiV 


02 


1 2 


1 0 








vj ▼ 


0 2 


i & 

lrU 


1 Q 








HAI 




1 


J.O 








* 


20 kHz 


A.I1ACUJ 






■ 




flu 


60 kHz 




\nxcQ) 






- 


MB* 


10 kHz 














60 kHz 


y jaw a v 

(fixed) 


(fixed) 








4B3 


70 kHz 


(fixed) 


(fixed) 


10 


2 


6 Koj 




0.5 




4.3 






Kio 


c 2 V 


0J 


1.0 


5.0 








qAi 


60kHz 


0.70 


1.4 






K„ 




20 kHz 


1.3 


4.5 






Km 


qsi 


40 kHz 


0.90 


2.9 








qaa 


80 kHz 


0.55 
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Error values are calculated from the scattered results of analysis applied to a series of 
simulated data, 

* ' * ,. * t "* 4 - . , h - ■ " i *. i 

r * •••• > i b ••. . j , - ' 

r _ _ ^ * " ■ ^ t V 

1 ■ « ' " ' '- L _^ J ' - " ^ , r 

In Table 4 the relative deviation of mean values (i.e. bias) of estimated parameters are 

'. '- ■ " ' ' ' *■-...'■'. It'., ■ !' ' i4 . ... ■ . ■ . .■ 

presented for 2D-FIDA (present invention) and 2D-MAFID (prior art); In each case bias is 
determined from analysis of a series of thirty simulated random distributions of count 
numbers; Three cases were analysed. In the first case, models used in data simulations and 
data analysis were identical. In the second case, the distributions of count numbers were 
simulated assuming that particles of the" second species are not equivalent but being 
distributed by their individual brightness with a relative half- width of 20 percent. This 
phenomenon was intentionally ignored in analysis, however. Of course, applying a slightly 
inadequate model for analysis produces bias of estimated parameters < The third case is similar 

, t * r * . • - | - . 

to the second one except the relative half-width of the individual brightness distribution of the 
second species is 50 percent, which is a usual value for vesicular preparations. It is worth 
noting that methodological deviations are noticeable when mapping weighted residuals of 2D- 
FBDA in cases two and three, but 2D-FEDA still returns meaningful results/ 2D-MAFID is a 
more sensitive method to model roughness. 



Error values arc calculated from the scattered result* of analysis applied lo a 
simulated data. 



In Table 4 the relative deviation of mean values (i.e. bias) of estimated parameters are 
presented for 2D-FIDA (present invention) and 2D-MAFID (prior art); In each case bias is 
determined from analysis of a series of thirty simulated random distributions of count 
numbers, Three cases were analysed, In the first case, models used in data simulations and 
data analysis were identical. In the second case, the distributions of count numbers were 
simulated assuming that particles of the second species arc not equivalent but being 
distributed by their individual brightness with a relative half-width of 20 percent. This 
phenomenon was intentionally ignored in analysis, however. Of couree f applying a slightly 
inadequate model for analysis produces bias of estimated parameters. The third case is similar 
to the second one except the relative half-width of the individual brightness distribution of the 
second species is 50 percent, which is a usual value for vesicular preparations. It is worth 
noting that methodological deviations are noticeable when mapping weighted residuals of 2D- 
FIDA in cases two and three, but 2D-FTDA still returns meaningful results. 2D-MAFID is a 
more sensitive method to model roughness. 
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Table 4 



Parameter 


True value 


Bias of 2D-FIDA, percent 


Bias of 2D-MAFTD, percent 


specification 


(selected) 






r 










* 


Case 1 


Case 2 


Case 3 


Case! 


Case 2 


Case 3 


c,V 


1.0 


+0.1 ±0.2 


-0.35 


-2.0 


+0.1 ±0.2 


-6.5 


-28 


c 2 V - ; V ; - 


0.05 


+0-6 ±0,4 


-2.0 


-10.6 


+1.8 ±0.6 


-15.6 


-69 


Qai 


20 kHz 


-0.1 ±0.2 


+0.8 


+5.0 


-0.3 ±0.3 


+10.0 


+58 




100 kHz 


-0.3 ±0.4 


-0.7 


-11-2 


■ * - ,* " 
i n ■ 

-0.9 ±0.3 


+5.5 


+41 


qai 


1kHz 


-1.2 ±0.5 


+2.6 . 


+20.6 


-6.7 ±3.0 


+75 : . 


+563 




200 kHz 


-0.3 ±0.4 


+1.2 


+1.3 


-0.9 ±0.3 


+1 1.6 


+99 



BRIEF DESCRIPTION OF FIGURES 

Figure 1 shows a count number distribution obtained for a solution of the dye 
tetramethylrhodamine (a) and residual curves corresponding to different fitting procedures 

Figure 2 shows simulated distributions of the number of counts for a case which models a 
binding reaction of a labeled ligand to vesicles. 

Figure 3 illustrates theoretical errors of the estimated parameters c and q of a solution of 
single species, depending on the value of q< 

Figure 4 illustrates theoretical errors of the estimated parameters c and q of a solution of 

. ■• ■ . ■ 

single species, depending on the value of c, 



5 illustrates theoretical errors of the estimated parameters of c and q of a solution of 
sing) species, depending on the value of T. 



f 
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Figure 6 illustrates theoretical enors of the estimated parameters c and q of a mixture of two 
species, depending on the value of T. 

Figure 7 illustrates theoretical errors of the estimated parameters c and q of a mixture of two 
species, depending on the ratio of q% to 

Figure 8 illustrates theoretical errors of the estimated parameters c and q of a mixture of two 
species, depending on concentrations. 

Figure 9 illustrates an apparatus adapted for use in performing the method according to the 
present invention. 

In Figure 10A, the calculated photon count number distributions, P(n), for cases having 
identical mean count number ft but differing by the composition of the sample are plotted. 

' * ■ ' ■ ■ ,.'.'*',. 
Figure 10 B illustrates the distributions of the number of photon counts of pure solutions of 

rhodamioe 6G (Rh6G) and tetramethylihodamine (TMR); as well as a mixture of these two 
dyes. 

' * < . 

The exemplary residuals for Rh6G are shown in Figure 10 C 



Figure 10 D shows the results of an TTR analysis (inverse uansformation with the help of 
linear reguiarization and constraining concentrations to non-negative values) applied to the 
curves of Figure 10 B. 



Figure 11 shows an APR analysis of hybridized (A-Q and restriction enzyme cleaved (D, E) 
labeled oligonucleotides. The curves result from a set of 20 individual 10 s measurements 
which show variations among each other of less than 10 %. 



Figure 12 illustrates hybridization and restriction enzyme cleavage of different combinations 
of labeled and unlabeled oligonucleotides. 



Figure 13 is a graphical presentation of a joint distribution of the number of photon counts 
, measured ! r a solution of TAMRA. 

I . . . .. • ■ V ■ . . - ' 

Figure 14 is a graphical presentation of weighted residuals of a joint distribution of count 
numbers obtained from a mixture of TAMRA (5 ? -(6^arboxytetramethylibodamine) and 
RRX (rhodamine red X). 

Figure IS illustrates joint distributions of the numbers of photon counts for the "parallel" and 
"perpendicular" polarization components of fluorescence measured for equal theophylline 
concentration of 2 nM, but different antibody concentrations. 

- . . ■ - ' l' - J ■ - * ^ * ♦ * - , . _ t r "... . ^ h . ., . ! . _ ^ * 

* ' ' ~ i * ' 1 ■ ' ' " ' * . " * ' ' * ' i " 1 

Figure 16 illustrates the results of 2D-FDDA applied to data from samples of equal 
theophylline concentration, but different antibody concentration. : - 

r. r * . . 
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Figure 17 illustrates the principle of a two-color 2D-FIDA experiment with multiple binding 
sites per vesicle. . 

Figure 18 illustrates joint distributions of the numbers of photon counts in "green" and "red" 
measured in conditions of high (A) and low (B) degree of binding of SMS- 14 (somatostatin- 
14) to SSTR-2 (human type-2 high-affinity somatostatin receptor). 

Figure 19 shows the competition curve of the binding reaction of SMS- 14 to SSTR-2. 

t ** * ■ i 

Figure 20 illustrates ah embodiment in which numbers of photon counts [ni] subject to 
determination of a distribution function P(n) in step b) are derived from numbers of photon 
counts in primary time intervals (Nj) by addition of numbers of photon counts from primary 
time intervals according to a predetermined rule. 

; ■ ' j ■ j ■ - - 

Figure 21 illustrates a further embodiment in which numbers of photon counts {n,} are 
derived from predetermined primary time intervals according to a rule in which primary time 
intervals are separated by a time delay. 



DETAILED DESCRIPTION OF FIGURES 

Reference is now nude to Figure 1 which ahow, a count number distribution obtained for a 
solution of the dye tctramethylrhodamine (a) and three residuals curves corresponding to the 
best fit obtained with Eqs (*), (3», and Q» below (b). The concentration of the aqueous 
tetnimethylrhodaraine solution was about 10 " M. Primary time intervals were 40 us. The data 

collection time was 60 s. 

The most widely used spatial profile model in FCS is the three-dimensional Gaussian profile 
with a single parameter of shape, the axial dimension ratio «n longitudinal and radial 
directions- Residuals curve with open circles illustrates the fit quality obtainable with the 
Gaussian profile. There are large systematic deviations in residuals. What is a sufficiently 
flexible model for fitting FCS data has turned out to be a rather inflexible and inadequate 
model for FID A. 

A model of the sample profile which has yielded a better fit of the measured distribution 
P{n) is Gaussian-Loientzian (open squares), but this model still lacks flexibility. According 
to Eq. (8), a certain function of the spatial brightness B is integrated over the volume. In other 
words, it is a relationship between B and V, characterizing a given spatial brightness profile in 
F1DA. For example, the Gaussian profile yields the relationship 



dV r 

04 VJt, 



dx * ' (36) 



-where x = -In B. The Gauss-Lorentzian profile yields the relationship 



(37) 



Both of the relationships are rather inflexible, i.e., they do not provide any spatial shape 
parameters to adjust the theoretically calculated distribution to fit the measured data. 

When looking for sufficiently flexible models to fit experimental data, it might be useful to 
apply the following relationship: 



4V 
tlx 



(38) 



There is a formal rather than a physical model behind Eq. (38); the fit quality obtainable with 
Eq. (38) is illustrated by the filled squares curve: 

Referring now to Figure 2, simulated distributions; of the number of counts are shown for a 
case which models a binding reaction of a labeled ligand to vesicles. Primary time intervals 
were 46 us, values of the spatial parameters of Eq. (38) a, = -0.4; flJ = 0.08V background count 
rate b= 1.0 kHz, data collection time 4 s. Curve if ligand" corresponds to a species of c = 6.0; 
q .•= 6.0 kHz/particle; oy= 0. Curve ..vesicles" corresponds to a species of c = 0.05; ? = 300.0 
kHi/particle; d, = 150.0. Curve „ mixture" corresponds to their mixture. Concentrations and 
specific brightness values have been selected to model a characteristic situation in drug 
screening. Fitting of curve (c) returns the values of the five parameters characterizing the 
given ..sample'- with statistical errors which are mostly between 3.5 and 6 percent, except the 
error of a, of vesicles which is 13 percent. If, however, a, of vesicles is fixed in fitting, all the 
statistical errors are below 4 percent 



For the fastest data simulation algorithm, one may calculate the expected distribution and 
generate a random Poisson number of events for each value of n indepenctently^^s a 
cosmetic error, the total number of events may slightly deviate from the pre-given 

number M. A slower but a straightforward data simulation algorithm is the generation of a 
random configuration of particles in volume elements contributing to fluorescence, the 
calculation of the classical light intensity corresponding to the given configuration of particles, 
and the generation of a random Poisson number corresponding to this intensity, as a simulated 
number of photon counts. This procedure is repeated M times to obtain a simulated count 
number distribution. 



* 



Reference is now made to Figure 3 which illustrates the theoretical errors of the estimated 
parameters c and q of a solution of single species, depending on the value of q. The following 
values of experimental parameters were selected: c = 1.0; T = 20|ts; a t « -0.4; ai = 0.08; b » 

1.0 kHz; data collection time 2 s. 

Generation of a random P ision number is the f Mowing. For a given expected value of events £". a simulated 
number of events 5 is determined from a routinely generated random number /f between 0 0 and 1.0 through 

inequations £ Poisson(i: £) < R < X PoissoniiiE).;., 
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Figure 4 illustrates the theoretical errors of the estimated parameters c and q of a solution of 
single species v depending on the value of c The following values of experimental parameters 
were selected: q = 60 kHz/particle; T = 20 jis; a x = -0.4; a 2 = 0.08; b = lo kHz; data 
collection time 2 s. 




Retelling now to Figure 5, theoretical errors of the estimated parameters c and q of a solution 
or single species are shown, depending on the value of T. The following values of 
experimental parameters were selected: c = 1.0; q ■ 30.0 kHz/particle (upper graphs): q = 60.0 
kHz/particle (lower graphs); a x - -0.4; a 2 = 0.08; b = 1 ,0 kHz; data collection time 2 s. 

Figure 6 illustrates the theoretical errors of the estimated parameters c and q of a mixture of 
two species, depending on the value of T. The following values of experimental parameters 
were selected: a - 0.1; c 2 = 2.0; = 200.0 kHz/particle; q 2 = 10.0 kHz/particle; a, = -0.4; a 2 
~ 0.08; b = 1.0 kHz; data collection time 10 s. Note that the optimal value of T for the 
detennination of the parameters of the brighter species is lower than that of the darker species. 

Reference is now made to Figure 7 which illustrates the theoretical errors of the estimated 
parameters c and q of a mixture of two species, depending on the ratio of q 2 to qi. The 
following values of experimental parameters were selected: q { = 50.0 kHz/particle; Ci ~ ci = 
0 ,5, 0i = -0.4; a 2 - 0.08; b = 1.0 kHz; data collection time 40 s. 

Figure 8 illustrates the theoretical errors of the estimated parameters c and q of a mixture of 
two species, depending on concentrations. The concentrations were changed synchronously, c\ 
- a. The following values of experimencal parameters were selected: q\ = 75.0 kHz/particle; 
q 2 = 25.0 kHz/particle; a { = -0.4; a% = 0.08; b = 1.0 kHz; data collection time 60 s. Note that 
an optimal concentration exists at about one particle per sample volume. This is generally 
true, except if less than three parameters are to be determined. 



Reference ts now made to Figure 9 which shows one embodiment of an apparatus adapted for 
use in performing the method according to the present invention. Apparatus 10 comprises a 
laser 12 which serves as a light source for illuminating the sample by a bundle of coherent, 
monochromatic excitation radiation 14. Excitation radiation 14 is paralleled by a lens 16 and 
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reaches a dichroie mirror 20, Preferably, the angle between the optical axes 18 and the 
dichroie mirror 20 is 45°. The dichroie mirror 20 reflects the excitation radiation 14 in 
direction of an objective lens 22 having its focus 24 within a sample volume .26. Sample 
volume 26 and objective lens 22. arc preferably separated from each other by a transparent 
cover glass 28i e. g* by the bottom of a commercially available micro^titer plate which houses 
the sample. The sample preferably includes fluorescently labelled molecules or other particles. 
Due to excitation by an appropriate excitation radiation 14, the molecules or other particles 
present in the sample emit radiation 30. Emission radiation 30 passes the objective lens 22 and 
r reaches the dichroie mirror 20 which is transparent for emission radiation 30. Thereafter, 
emission radiation passes a filter 34 and a collimator lens 36 on the optical axes 32 A pinhole 
38 is situated in the focus of collimator lens 36: Emission radiation 30 passing the pinhole 38 
reaches a farther lens 40 and, thereafter, is detected by the photo-detector 42. Within the 
pathway of emission radiation 30, in particular between dichroie mirror 20 and photo-detector 
42, an opaque means 44 is provided through which a central part of the emission radiation 30 
cannot pass. This central part of the emission radiation 30 stems from areas on the optical 
axes 32 in front of or behind the focus 24 of the excitation radiation 14. Only emission 
radiation 30 that stems from the focus 24 or its direct neighbourhood passes the pinhole 38 
and reaches photo-detector 42. Instead of placing an opaque means 44 within the pathway of 
emission radiation 30, the pathway of excitation radiation 14 is also suitable for positioning an 
opaque means 44. In particular, an opaque means 44 can be positioned between laser 12 and 
dichroie mirror 20. Use of an opaque means 44 as described in detail herein improves the 
signal-to-noise ratio. 

" ' . * ' * ' ' • 1 l^ : ' " 1 

In Figure 10A, the calculated photon count number distributions v for five cases of equal 
mean, count rate, #T« L0 are shown. The open symbols correspond to solutions of single 
species, but with different values of the mean count number per particle qT . The solid line is 
calculated for a mixture of two species, one with = 0-5 and the other with qT = 8,0 . 
Obviously, the curves differ from each other considerably. The important point is that FIDA 
according to the present invention can unambiguously separate the contributions of the 
individual species. 

. . r * ■ . - . .. .. 

Reference, is now made to Figure 10 B which illustrates the distributions of the number of 

• • ■ . ' - ■■ .. - . ' ' T . " * -» • • 

photon counts of pure solutions of two different dyes, 0.5 nM rhodamine 6G (Rh6G) and 1.5 
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oM letrametbylrhodamine (TMR), as well as a mixture of the two '(0.8 nM TMR. 0.1 nM 
Rh6G). The main equipment is a confocal microscope (ConfoCor®; EVOTEC BipSystems 
and Carl Zeiss. Germany) routinely used for fluorescence correlation studies. An attenuated 
(to about 800 nW) beam from an argon ion laser, wavelength 514.5 nm, is focussed to a spot 
of approximately 0.5 jim radius, which is twice the size of a spot in usual FCS experiments, 
and results in a diffusion time of approximately 200 us/for rhodamine 6G. the excitation 
intensity has generally been kept lower than or equal to a level characterized by about 15 
^nt amplitude of 4e triplet s 

detected through a pinhole on the focal plane of the microscope using an avalanche 
photodiode detector SPCM-AQ 131 (EG&O). The distributions of the number of photon 
counts were measured at T = 40 us dwell time. Data collection lime was 50 s. These 
distributions serve as input data for FEDA. The results of a multisroinponent fit analysis are 
given in Table 5 below with % l values calculated according: 

/ _Ew.[p(«)-/»wr (39) 

■ ftp ~ n jit , 

where n p is the length of the measured distribution > and n p is the number of fit 
parameters, In said mulii-compbneot analysisooo fits the measured distribution function of 
number of photon counts, assuming a certain number of fluorescent species, and estimates 
unknown concentration and specific brightness values. 



V. ' I ■! 

Sample 








gVMi) 




RH6C 


-0.380 ±0.009 


0.077 ±0.003 


0461 ±0.003 


107.2 ±6:8 : 


0.97 


TMR 


-6.427 ±0.032 


0.084 ±0014 


1.517 ±0.012 . k 


36^6 ±0 29 


0.81 


Mixture 


-0 380 (rued) 


0.077 (fiwd) 


0,103±0.012 :• , 


109i±4.0 

• - 


0.77 
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0738 ±0.011- ^ 


37.4 ±1.0 ; 


. 1 . 



In all cases, the background count rate was fixed to, the value of t ,05 kHz, as measured with 



de-ionized water, ai and ay are pre-normaiization values when a \ is fixed to 1-0. The 
theoretical statistical errors in Table 3 correspond to theoretical weights; 

1 1 - ' + : , - - 

" i ' - . , -/ ' ••• . ' -■ "7 • • . ', '- 

.' -•-|V B =^ -J-.-. (40) 

This formula is derived under a simple assumption of TV independent measurements of the 
number of photon counts; to reality, consecutive measurements are correlated, therefore the 
errors of estimated parameters returned by the fitting algorithm underestiinate real statistical 
errors We have empirically determined from a separate series of 30 to 200 measurements that 
these statistical errors are greater than theoretical ones by a factor of about three. 

The exemplary residuals for Rh6G are shown in Figure 10 C. 

The results of ITR analysis are shown in Figure 10 D expressing the distribution of the 
number of particles as a function of their specific brightness. Such ITR analysis is an inverse 
transformation realized with the help of linear regularization and constraining concentrations 
to n on- negative values. ITR is a valuable tool especially for samples from which either no a 
priori information about the sample composition is given cur where the sample composition is 
heterogeneous. The dashed lines correspond to the solutions of the single dyes (Rh6G and 
TMR) and the solid line to their mixture. The ordinates give the mean number of particles 
within the cpnfocal Volume element Oeft: single dyes; right: mixture). For equipment details 
see description of Figure 10B. The ideal outcome would be single J -peaks for the solutions of 



single species and two tf-peaks for the mixture. In reality, the width of the ITR output spectral 
peaks is determined not only by the true width of the distribution of specific brightness values 
but also by the accuracy of input data and the particular realization of the linear regularization: 
simply, if a broad spectrum fits experimental data as well as a narrow one then ITR prefers the 
broad one. 




As a further example, the method according to the present invention has been applied to study 
ihe hybridization of 5-(6-carboxyieLramethylrhodarnine (TAMRA)Habeled 40mers with 
either labeled or non^labeled complementary oligonucleotides and the subsequent symmetrical 
cleavage of the DNA hybrid by the restriction endonucleases Hind /// and Kpn I. The specific 
oligonucleotides used in this study were 

TAMRA-AAGAAGGGGTACCTTTGGATAAAAGACAAGCT^ 
Oligo A) and 

TAMRA-ACGGGAAAAGCTTCTCTTTTATCCAAA '-TAMRA- 
Oligo B). 

They were purchased in HPLC pure quality from Applied Biosystems (Weiterstadl, Germany). 
All measurements were carried out on the above described FCS reader at excitation/emission 
wavelengths of 543/580 nm using a 4 mW helium/neon laser (Uniphase) attenuated to 
approximately 300 mW. The water background was found to be below 800 Hz. For the 
measurements, sample aliquots were diluted to 1 nM and 20 \i\ were assayed in a 8- well 
chambered coverglass (Nalge Nunc) at room temperature. The hybridization reaction was 
performed in 70% formamide containing lOoiM TnsfHCl buffer (pH 8,0), ImM EDTA, 
0.2mM NaCl and an oligonucleotide concentration of 0.5^M. Denaturacion was at 95°C for 2 
min and subsequent hybridization was at 55-60°C for 40 min in accordance with the 
optimized temperature Tx which is 10-15 degrees below the melting point Tm (Heating block, 
Techne). Tm was calculated as follows: Tm = 81.5 +16.6 (logl0[Na+]) + 0.41 (%G-+C) - 
600/N (N = 40; %G+C = 42,5), Restriction digest analyses of the hybrid DNA were performed 
by the restriction enzymes Hind III and Kpn I The restriction site was chosen in order to 
obtain fragments of different size. The cleavage reactions were performed in 3.3mM 
iris/acetate (pH 7.9), 1 mM magnesium acetate, 6.6. mM potassium acetate and O.lmg/ml 
bovine serum albumin at 37 °C for lh. The reaction course is characterized by significant 
shifts of the fluorescence intensity per molecule within a range of one order of magnitude 
(Figures 11 A • E). It is important to note that the resulting double product peak corresponds 



to one single hybridization species in which no none-hybridized species are delectable. This 
could be easily demonstrated in hybridization/cleavage experiments using labeled/nonlabeled 
oligonucleotide combinations where signals with different brightness compared with the 
double labeled hybrid have been obtained (Fig. 12 A f B). One explanation for such a double- 
peak structure could be conformational fluctuations resulting from transitions between two 
conformational states. This behavior is described using a three-state model of the 
conformational dynamics vyith a polar* a nonpolar, and a quenching environment of the label 
(Eggeling et al. p Proc. Natl. Acad. Sci- USA 95: 1556 - 1561). From, the present data it can be 
concluded that the cbnfoirnatibnal fluctuations are significantly dependent on the DNA length 
and the number of labels. The small digested DNA fragments exhibit a single brightness and 
the concentration ratio of the educt double peaks shifts toward one component when one-label 
hybrids are studied. As a further result, small fragments show low molecular brightness and 
large fragments correspond to high brightness values (Q2 *■ Q3» Qi Q*). , 

Based on symmetrical cleavage sites for both enzymes DNA fragments of equal size and 
subsequently very similar molecular brightness are to be expected. Figures 12 A and B 
clearly illustrate that such a result has been observed and the theoretical relationship Qi + 02 = 
Q3 + Q4 is experimentally confirmed. Restriction digest analyses of the hybrid DNA were 
performed with the restriction enzymes Hind 01 and Kpn 1 either alone or in combination. The 
restriction sites lead 10 symmetrically cleaved fragments (2x12 mcrs, 1x16 mer for double 
digest). [•] corresponds to labeled DNA, Qdh»: intensities of DNA hybrids, Q;: intensities of 
cleavage products. The bulk of all the measurements of cleavage products using single and 
doubly labeled DNA structures makes it possible to classify all observed intensities with 
individual structures in a single measurement - a property of the method of the present 
invention which is not possible using any other known analytical method. 
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Reference is now made to Figures 13 and 14. As the central optical part of the equipment for 
2D-FIDA* a confocal microscope is used, like in FCS spectrometers (Koppel et al., Biophys. 
J. 16: 1315 * 1329, 1976). For excitation of fluorescence, a beam from an Ar or green He-Ne 
laser is attenuated by neutral Biters, passes a beam expander and is directed to the microscope 
objective by a dichroic mirror. In a number of experiments with slowly diffusing particles, 
beam scanning in combination with sample scanning is used, as a tool known from laser 




scanning microscopy. Fluorescence is collected by the same objective through the dichroic 
mirror, and is focussed to a confocal piahcle which serves to reject the out- of- foe us light 
Resolution in the longitudinal direction is additionally improved by using a concentrical 
opaque spot closing about a quarter of the aperture in the fluorescence collection path of the 
microscope (see description of Figure 9). The light which passes the pinhole is divided by a 
beamsplitter for detection by two detector*. Depending on the general type of a 2D-FIDA 
experiment, the beamsplitter is either a polarization cube, or a dichroic mirror, In the fust 
case, a common spectral band-pass filter is used, while in the case of two-colour FIDA, each 
detector has a different band-pass filter in front of it The photon counting detectors are silicon 
avalanche photodiode modules SPCM-AQ-131, EG&G Optoelectronics, Canada. The TTL 
pulses from the detectors are counted by a two-channel counter, constructed by EVOTEC 
BioSystems AG (Hamburg, Germany) as a plug-in card of a computer. The count number 
distributions are calculated during reading data from the 32 MB onboard buffer, which itself is 
an online process. By feeding the detector outputs to a correlator, FCS measurements can be 
performed in parallel with FEDA experiments. 




The levels of background count rate for both detectors are determined by a separate 
experiment on bi-distilled water. The main contributor to the non-fluctuating background light 
intensity is Raman scattering from water. 



The radius of the monitored sample volume can be adjusted by selecting an appropriate 
expansion factor of the original laser beam. The focal beam radius of about 0.6 fun is used 
yielding diffusion times for simple organic dye molecules (e.g., TAMRA) of about 260 jxs, 
which is considered long compared to the 40 us dwell time of counters, so that the 
approximation of constant molecular brightness during the counting interval is valid The 



excitation intensity is adjusted as a compromise between a high count rate per molecule and 
low population of the triplet state. The triplet state population was kept at about 15 percent 
level. Higher triplet population values might significantly distort the apparent spatial 
brightness profile. The population of the triplet state was measured with the correlator (ALV- 
5000, AJLV, Langen, Germany) during the set-up of a series of experiments. 

For methodological test experiments, two different dyes were selected, carboxytetramethyl- 
rhodamine (TAMRA) and rhodamine red X (RRX), These dyes have different emission 
spectra, as well as different extinction coefficients at the laser wavelength of 543-5 am. In 
these experiments, a wideband 40/60 beamsplitter was used in front of the detectors. The 
spectral filter of the ,jed" channel has the center wavelength of 605 am and FWHM of 50 nm 

while the corresponding figures for the Pf yellow- green" channel are 575 nm and 30 am. The 
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dyes were diluted in distilled water so that the average number of molecules in the observation 
volume was in the range of 0.5 to 2.0 t corresponding to concentrations between 0.23 and 0.92 

. nM. 

For each experiment, about 20 |il of the sample solution was placed on a coverslip separating 

■ J * ' *. ' - ' \ ■ - - „ > ■ j 
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the sample from the water immersion objective (Zeiss C-Appchromat 40x1.2 W Konr), Each 
distribution was collected for 60 seconds. Each of the dyes was measured separately, but also 
mixtures of the dyes with different concentration ratios were measured. The parameters 
describing the spatial brightness profile were determined from an experiment on TAMRA and 
were fixed in subsequent analysis of other samples at values a x = -0405 and a, =0.0772. 

' ' *' . ■ ' f ■)■. . -r • . '< 

Mixtures were measured and analyzed in order to see how well the method returns values of 

. • ■ ■ ■ * ' ■, ■ • * . . ' ■ . 

specific parameters of the two species. 
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As an example/ Figure 13 visualizes a count number distribution measured for a pure 
TAMRA solution. Values corresponding to ihe *-axis are logarithm values of the probability 
to obtain a given pair of count numbers m and /12. Fitting of the distribution returns mean 
number of panicles. rV = 1.139 ±0.005, specific brightness for the „red" channel qi = 79.4 
±0.5 kHz, and for the ^yellow-green" channel q% = 50.9 +0-3 kHz, In Table 6, results of 
analysis of test experiments are presented. Samples have not been specified by concentration 
values calculated from dilution factors of the preparation because adsorption of dye molecules 
to glass surfaces influences real concentrations significantly. Because of the same reason, 
concentration Values are not well reproduced from sample to sample; a shift of concentration 
values from realization to realization is also sometimes observable. Specific brightness values 
are well reproduced, however. Values of statistical errors presented are theoretical values 
corresponding to the assumption of ^correlated measurements, multiplied by three, which is 
an empirical factor. In addition to statistical errors, deviations from sample to sample of a 
modest size are also noticeable. 



• 



Table 6 



Sample 


Realizaiio 


Mean number of 


Specific brightness 


Specific brightness 




d number 


molecules per 
sample volume, cV 


in „red", qj, kHz 

• i" ' 


in,, yellow-green", 
q2,,kHi 


TAMRA 


1 


1.128 ±0.005 


79.6 ±0.4 


51.3 ±0.3 




2 


1.139 


79.4 


50.9 




3 


1.160 


79.4 


508 






1.171 


79.0 


50.6 


RAX 


I 


1.892 ±0.009 


48.0 ±0.3 


15.63 ±0.09 




2 


1.859 


48.0 


15.63 




3 


1.850 


47.3 


15.44 




4 


1.811 


47.9 


15.58 





-;3 J 0 i -j . i ; 



- -6-99- V.'r 



TAMRA and 
RRX 



0.99 ±0.05 
1.14 ±0.05 
1.01 
1,08 . : 
0.99 . 
1 .03 
0.97 
1.03 



d^m/mmmM^ 22i 
II 

36 ■ , T : 

■■ j- + 

.'V ■ - 

• ' . . < ' ' ' . ■, 

78.0 ±1.0 
^ 49.1 ±1.1 
76.8 
.. 49.5 

', * - ** -V " 

- ■ .f ■ I . 

77.9 

493 

■ y 76:9. . 
50.0 
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+49,89 
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51.8 ±1.1 

15.0 ±0.9 

51.1 

15.1 

51.8 

15.1 

51.7 1 

' T - 

15.5 



TAMRA and 


- I .• 


0.93 ±0.07 


74,6 ±1.3 ,,' 


49.6 ±1.5 


RRX, new 


» - ' - ■ 


\y . 2.51 ±0.07 " 


■;: ■ 47.6 ±o.8 


14.7 ±as 


sample 




0.98 




49.2 ■• '-' V;- ■ 






2.38 


46.8 / 


14.2 


H ' . . - i - 


- 3 


0.96 

-I l 


75.0 ' 


• 49.4 


* . * ■ J" 
1 

«■ * 




. 2.25 . 


48.0 


" ' 14.7 




i: 4 - 


0.94 


..../, 76 a . ,. 


■ 4*9 . 




t 


2.16 .•; - 


;.V-. 47.6. ■ 'TV 


14.8 • 




Ic is worth noting that distributions measured with mixtures are qualitatively different from 
those measured for pure dyes. Figure 14 visualizes weighted residuals of fitting a distribution 
measured with a mixture of TAMRA and RRX The upper graph corresponds to the adequate 
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analysis when two species were assumed to be present; residuals are scattered quite randomly 
and uniformly. The lower graph corresponds to the assumption that only a single species is 
present; there is a significant difference between the measured and the calculated distribution 
in this case, m is the count number obtained by the „red" detector, and n% is the count number 
obtained by the „yellow-green* 4 detector. 



€ 



In traditional fluorescence polarization studies, average intensities of two polarization 
components of fluorescence are direcdy measured. Fluctuations of the intensities are not of 



direct interest there but are considered rather as a source of statistical errors. Changes in the 
fluorescence polarization values of a sample containing a fluorescently labeled binding partner 
reflect changes in molecular volume and, hence, provide direct information on equilibrium 
binding. Fluorescence polarization measurements can also be performed in real-time, allowing 
the kinetic analysis of association and dissociation reactions. One of the most widely-used 
fluorescence polarization applications is the competitive immunoassay used for the detection 
of therapeutic and illicit drugs. The method of fluorescence polarization has been used for 
clinical immunoassays for more than a decade. The homogeneous FPIA (fluorescence 
polarization in immunoassays) has well-accepted advantages over conventional heterogeneous 
immunoassays like R1A or ELBA. However, it fails if multi-binding step reactions are to be 
investigated because the separation of individually polarized species is impossible. Therefore, 
ligand binding curves only demonstrate the overall decrease of polarization, meaning the 
mechanistic binding constants cannot be determined. Further limitations are seen in sample 
volume as well as in mass restrictions. 



With a two-dimensional fluorescence intensity distribution analysis according to the present 
^ invention, the full content of information usually buried in fluorescence anisotropy can be 
utilized, thereby overcoming the limitations mentioned above. 2D-FDDA directly determines 
two specific quantities per each fluorescence species in one measurement: the fluorescence 
intensity per molecule and the anisotropy of a given model. Based on this supplementary 
information, the delineation of all participating species and even the quantification of the 
binding behavior is possible. 2D-FTDA anisotropy is an ideal tool for the quantitative 
description of systems exhibiting multiple binding steps, aggregation and multimerization 
phenomena, 
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Theophylline therapy has been a cornerstone of asthma therapy for several years and, 
therefore, there is a strong demand for assaying and fine-tuning the theophylline level in 



serum. To demonstrate the capabilities of the present invention, the binding of theophylline 
antigens to ariu theophylline anUbodies ^ a polyclonal anti-thcophylline serum has been 
investigated. The antigens were labeled with 5%cafboxytetraniethylrhodamine (TAMRA) with 
and without an incorporated spacer. In classical FP analysis these conjugates exhibited low 
inisotrbpy values (0.037 and 0 055 respectively) upon interacting with the antibody. Therefore 
they form a critical basis for illustrating the sensitivity of the present invention. 



For the binding experiments the stock solutions of antibody and antigens were diluted in a 
PBS buffer with 6.05% Tween,20. After rnixing the compounds, the mixture was incubated 

■ . - - • •." ' - ' ,- ' , . . 

for 30 minutes at room temperature. Matching the spectral properties of the conjugates the 
system was excited at the 543 run line of a He/Ne laser (Uniphase). As two examples, 
measured joint distributions of photon count numbers and results of 2D-FIDA applied to 
samples at different antibody dilution values but a constant ligand concentration of [L] =2nM 
arc illustrated by Figures 15 and 16. Andbody concentrations arc in arbitrary units referring to 
effective dilutions. Water background was below 1 kHz in each detection channel. 



In the two-color fluorescence intensity distribution analysis according to the present invention; 
two detectors are spectrally tuned to monitor fluorescence from two labels of different color, 

' . . ' ■■■»* * ' "'"-■.■„-' ' 1 j - ;• • ■ ' -.'.■.,'» ■ 

In the assay type described below, Ugand molecules are labeled in "green". Since each vesicle 
carries a high number of receptors, vesicles in samples with a low binding degree can be 
distinguished from vesicles in samples with high binding degree by a significantly higher 
specific brightness in "green rt . Vesicles are additionally stained in "red". Specific brightness 
of vesicles in "red!' is not altered by binding of ligand molecules, but staining in "red" is a 



means to increase contrast between free ligand molecules (which arc nearly invisible in •■red") 
and vesicles (which in the case of extremely low binding may be of nearly the same brightness 
in "green" as free ligand molecules). Contributions from the two fluorescent species of a 
single sample to the measured 2-dimcnsional distribution of the numbers of photon counts are 
very different in this assay type indeed, and therefore the analysis is highly reliable. 



To demonstrate the advantages of a two-dimensional fluorescence intensity distribution 
analysis according to the present invention, the binding of TAMRA-labeled somatostatin-14 
£ (SMS 14-5TAMRA) to the human type-2 high affinity somatostatin receptor SSTR-2 (P. 
Schoeffter et al„ Eur. i. Pharm., 289, 163 - 173, 1995) has been chosen as biological system. 
The receptor was expressed by CCL39hsst2 cells. Vesicles were prepared from this cell line 
and stained by the lipophilic tracer DiC»(5). The binding reaction was carried out in 10 mM 
HEPES (pH 7.6). 5 mM MgCl,. 0.01% (w/v) fluorosurfactam FC-135, 1.33% DMSO in the 
presence of protease inhibitors. A schematical drawing of the principle of the assay is shown 
in Figure 17. 



For excitation of fluorescence of the two spectrally distinct labels, the 532 nm lice of a 
Neodyra:Yag laser attenuated to 250 nW and the 632 nm line of a He-Ne laser attenuated to 
25 jiW were simultaneously used. In order to minimize misalignment of the two laser beams, 
they both passed through a single optical fiber before being focussed by the microscope. Also, 
the collected fluorescence beam passed a single pinhole before being spliued for the two 
detectors. An optical band-pass filter with the center wavelength 590 nm, FWHM 60 nm, and 
another one with the center wavelength 690 nm, FWHM 40 nm were used in front of the two 
detectors, monitoring fluorescence from the two labels separately. Because vesicles are slowly 
diffusing particles, in this experiment an area of 0,05 mm 2 of each sample was scanned, using 



sinusoidal beam scanning of 25 tiz frequency, 100 fim amplitude in one direction, and sample 
scanning of 500 Jim per 8 s data collection time in the other direction. 

In Figure 18 two typical examples of count number distributions corresponding to a high and 
a low degree of binding are compared. In Figure 18 B. the vesicle bound ligand (SMS 14- 
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STAMRA) was competed off by the addition of a large excess of non fluorescent competitor 
(SR0M4, SIGMA; 1 \iM competitor, 3 nM total ligand). It is clearly visible that the 
distribution obtained in the absence of competitor (Figure 18 A) is more expanded in the ni 
direction than in the presence of competitor. The tlx direction represents the green 
fluorescence of the labeled ligand and, hence, is a measure for the Ligand binding. 

2D-FIDA according to the present invention was applied to each of the measured 
distributions, Using a two-component fit algorithm, both the concentration (particle number in 
the confocal volume) and the fluorecence brightness values for the free ligand and the ligand 
bound to vesicles were obtained. The vesicles are well separated from the free ligand even in 
the presence of competitor (low binding). Addition of competitor leads to a decrease in the 
green fluorescence of the vesicles and an increase in the concentration of free ligand. The fact 
that the ligand does not completely dissociate from the vesicles is due to unspecific binding 
which is dependent on the amount of vesicles used in the assay. 

Next, the competitor has been titrated from 1 pM down to 10 pM (10 measurements each) in 
order to record a dose-response curve and to determine the ECso value. As shown in Figure 
19* 2D-FIDA analysis leads to a typical sigmoidal competition curve with very low standard 
deviations. The calculated EC* value is 0.87 nM which is in good agreement with EC* values 
obtained using other evaluation methods (data hot shown). Thus, it is n t nly possible to 



differentiate between high and low binding but also small changes in the binding degree can 
be resolved with 2D-FIDA analysis according to the present invention. This is in particular 
important if one aims at the identification of assay inhibitors ("hits") in high throughput 
screening. In summary, 2D-F1DA anaJysis is a well-suited and highly reliable method for 
vesicle-based binding assays and has already been successfully applied in high throughput 
screening. 
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Reference is now made to Figure 20. According to a preferred embodiment, numbers of 
photon counts | n,} subject to determination of a distribution function P(n) in step b) are 
derived from numbers of photon counts in primary time intervals |Nj} by addition of numbers 
of photon counts from primary time intervals according to a predetermined rule. One might 
c.g. be interested in choosing numbers of photon counts {*] subject to determination of a 
distribution function P(n) which arc calculated from the numbers of photon counts in 

primary time intervals { Nj} according to the rule n, = £ N -(+ , , where M is an integer number 

expressing how many times the time interval in which {n,| is determined is longer than the 
primary time interval. Line N shows the primary time interval windows, Line a indicates the 
primary time windows chosen to calculate a, according to the rule. 

Figure 21 shows a further embodiment in which numbers of photon counts |nj are derived 
from predetermined primary time intervals according to a rule in which primary time intervals 
are separated by a time delay. In particular, the Mowing rule can be applied: 

M 

«* -EiN^+N^^), where M and L are positive integer numbers, [tk x \ are numbers 



of photon counts subject to determination of a distribution function P(n), and [Nj| are the 
numbers of photon counts in primary time intervals. 



CLAIMS 
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l A method for characterizing fluorescent molecules or other particles in samples comprising 

the steps of: • ; ' \ 

JL ■ . . *■ ...... ? . ■ 

a) monitoring fluctuating intensity of fluorescence emitted by the molecules or other particles 

in at least one measurement volume of a non-uniform spatial brightness profile by 
measuring numbers of photon counts in primary time intervals by a single or more photon 

detectors, • ; .. ' . "• "• 

b) determining at least one distribution function of numbers of photon counts, P(n) t from the 

measured numbers of photon counts, 

c) determining physical quantities characteristic to said particles by fitting the experimentally 
determined distribution function of numbers of photon counts, 

wherein the fitting procedure involves calculation of a theoretical distribution function of the 
number of photon counts P{n) through its generating function, defined as 

G(f) = £l" 

* _ ■* ■ J .-*'-." ^ " ' >' . ' ♦ 1 " * * * 
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2. A method according to claim 1 wherein the primary time intervals are consecutive intervals 
of equal width, 

3. A method according to one of claims l and 2 wherein in step b) numbers of photon counts 

fn; ] subject to determination of a distribution function P(n) are derived from numbers of 

photon counts in primary time intervals (Nj I by addition of numbers of photon counts from 
primary time intervals according to a predetermined rule* 

4. A method according to claim 3 wherein in step b) numbers of photon counts {n,} subject to 
determination of a distribution function P(ti) are calculated from the numbers of photon 

counts in primary time intervals { Nj) according to the rule n, = j^Nj^ , where M is an 

integer number expressing how many times the time interval in which (n, } is determined is 
longer than the primary time interval, 



51 A method according to claim 3 wherein in step b) numbers of ph ton counts (ii|) are 
calculated from predetermined primary time intervals according to a rule in which primary 
time intervals are separated by a time delay. 

6. A method according to claim 5 wherein in step b) numbers of photon counts (n;) subject to 
determination of a distribution function P(n) ait calculated from the numbers of photon 

counts in primary time intervals (NjJ according to the rule ^ B £'(N|^+N M((W ) ( 




where M and L are positive integer numbers. 

7. A method according to one of the claims 3 to 6 wherein in step b) a set of distribution 
functions P(n) is determined according to a set of different rules, said set of distribution 
functions being fitted jointly in step c). 




8. A method according to claim 7 wherein in step c) a set of distribution functions with 
different values of M and/or L are fitted jointly. 

9. A method according to one of the claims 1 to 8 wherein at least one of the physical 
quantities of step c) is concentration of particles* 

10. A method according to one of the claims 1 to 9 wherein at least one of the physical 
quantities of step c) is specific brightness of particles. 



1 1. A method according to one of the claims 1 to 10, wherein at least one of the physical 
quantities of step c) is diffusion coefficient. 



12. A method according to one of the claims 1 to 11 wherein the generating function is 
calculated using the expression G(£) = cxp[jdqc{q)jd y w{e { ^ m ^ -l)] f where c{q) is 

the density of panicles with specific brightness q % T is the length of the counting interval, 
and B(r) is the spatial brightness profile as a function of coordinates. 
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16. A method acooitUBg to any of the claims 1 to 15 wl 
is used for monitoring the mtensity offluore^^ 

, th* method according to any of the claims 1 to 

17. Use of a confocal apparatus for performing the tnetnoo bccoiwuib 

a i> f«r nmvicfinff excitation radiation (14), 

aradiationsaira 
an objecu™ ^ 

(26). . w - 

» detector (42) for detecting en^on radiation fc ™ 

volttme (26). and . . ^.^ r»\ ar 

excitation radiation O*)- 



the measurement 
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ABSTRACT 



A method for characterizing fluorescent molecules or other particles in samples comprising 
the steps of: 




a) monitoring fluctuating intensity of fluorescence emitted by the molecules or other particles 
in at least one measurement volume of a nonuniform spatial brightness profile by 
measuring numbers of photon counts in primary time intervals by a single or moit photon 
detectors, 

b) determining at least one distribution function of numbers of photon counts, P(n), from the 
measured numbers of photon counts, 

c) determining physical quantities characteristic to said panicles by fitting the experimentally 
determined distribution function of numbers of photon counts, 

wherein the fitting procedure involves calculation of a theoretical distribution function of the 
number of photon counts P(n) through its generating function, defined as 
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